Revising the predictions of inflation for the cosmic microwave background anisotropics 
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We point out that if quantum field renormalization is taken into account, and the counterterms 
are evaluated at the Hubble-radius crossing time or few e-foldings after it, the predictions of slow- 
roll inflation for both the scalar and tensorial power spectrum change significantly. This leads to a 
change in the consistency condition that relates the tensor-to-scalar amplitude ratio with spectral 
indices. A reexamination of the potentials </> 2 ,0 4 , shows that both are compatible with five-year 
WMAP data. Only when the counterterms are evaluated at much larger times beyond the end of 
inflation one recovers the standard predictions. The alternative predictions presented here may soon 
come within the range of measurement of near-future experiments. 



(N 

u: 
43 : 

Oh- 
i < 

o ■ 
> ' 

on : 

m . 
■ 

p ; 

o : 

ON . 
O ■ 



PACS numbers: 98.80.Cq 



A sufficiently long period of accelerated expansion in 
the very early universe is able to solve the questions 
raised by the standard big bang cosmology The hot 
big bang cosmology is an extremely successful theory. 
It explains the existence of the cosmic microwave back- 
ground (CMB) and its thermal nature, the observed ex- 
pansion of the universe, the abundance of light elements 
and the astrophysical fits for the age of the universe. 
However, it leaves without answer why our universe ap- 
pears so homogeneous and nearly flat at large scales. In- 
flation offers a natural answer to these questions and, at 
the same time, provides a predictive mechanism to ac- 
count for the small observed inhomogeneities 0] respon- 
sible for the structure formation in the universe and the 
anisotropics present in the cosmic microwave background 
(CMB), as first detected by the COBE satellite and fur- 
ther analyzed by the Wilkinson Microwave Anisotropy 
Probe (WMAP) satellite [3'] . Inflation predicts produc- 
tion of primordial density perturbations and relic grav- 
itational waves as amplifications of vacuum fluctuations 
together with a quantum-to-classical transition at the 
scale of Hubble sphere crossing. Primordial perturba- 
tions leave an imprint in the CMB anisotropies, which 
are, therefore, of major importance for understanding 
our universe and its origin. The potential-energy den- 
sity of a scalar (inflaton) field is assumed to cause the 
inflationary expansion, and the amplification of its quan- 
tum fluctuations and those of the metric are inevitable 
consequences in an expanding universe The metric 
fluctuations provide the initial conditions for the acoustic 
oscillations of the plasma at the onset of the subsequent 
radiation-dominated epoch. The detection of the effects 



of primordial gravitational waves in future high-precision 
measurements of the CMB anisotropies, as for instance in 
the PLANCK satellite mission [5(, will serve as a highly 
non-trivial test for inflation. Therefore, it is particularly 
important to scrutinize, from all points of view, the stan- 
dard predictions of inflation (as summarized for instance, 
in [6[) to be tested empirically. This is the aim of this pa- 
per. We point out that if quantum field renormalization 
is taken into account, as in the experimentally tested 
Casimir effect, the quantitative predictions of inflation 
change significantly, and may be tested in forthcoming 
CMB measurements. 

The scalar perturbations, which constitute the "seeds" 
for structure formation, are characterized by the power 
spectrum 
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where 4> represents the inflaton scalar field, which domi- 
nates the energy density during inflation. Here H stands 
for the Hubble rate H = a/a, where a(t) is the ex- 
pansion factor and dot means derivative with respect 
to the comoving time. The above expression is eval- 
uated at the Hubble radius crossing time tk (usually 
called "horizon crossing" time), where k/a(tk) = H. In 
the typical slow-roll inflationary scenario the homoge- 
neous part of the inflaton field <f>o(i) rolls slowly down 
its potential V((j>) towards a minimum. Both </> and 

H = J ^j^V(<po) are changing very gradually and this 

change is parameterized by the slow-roll parameters e, rj, 
where e = —H/H 2 , and r\ — e = <fr/(H<j>). These pa- 
rameters can be related to the derivatives of the inflaton 
potential e = {M 2 /2){V /V) 2 , rj = M 2 (V"/V), where 
Mp = l/y/8nG is the reduced Planck mass in natural 
units h = 1 = c. In the slow-roll approximation, e <C 1 
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and 1 77 1 <C 1, the scalar power spectrum turns out to be 



Pn(k) = 



2M 2 e(t k ) 



H{t k ) 
2ir 



(2) 



In addition, the power spectrum of tensor fluctuations is 
given by 



Pt{k) 



M 2 



( H(t k ) 
\ 2ir 



(3) 



and the tensor-to-scalar ratio is then r = Pt/Pji = 
16e. The power spectra are not exactly scale invari- 
ant, they can vary with k and this dependence is pa- 
rameterized by the scalar and tensorial spectral indices 
n s — 1 = dlnPiz/dhik, nt = dhiPt/dhik. Since these 
indices are related to the slow-roll parameters n s — 1 = 
— 6e + 2r] , n t = — 2e, one generates immediately the con- 
sistency relation r — — 8nt, which should be verified by 
any single-field slow-roll inflationary model, irrespective 
of the particular form of the potential. 

In this paper we shall reexamine, on the basis of gen- 
eral principles of quantum field theory in an expanding 
background 0, [I[ , the fundamental expressions (J2IS1) for 
the scalar and tensorial power spectrum. In doing this 
we shall also be led to modify the expressions for the 
spectral indices in terms of the slow-roll parameters and, 
therefore, to generate a new consistency relation. We 
will have all the necessary ingredients to reexamine the 
observational predictions of inflationary models and here 
we shall do that for some of the most significant models. 

The second factor in the fundamental relation ([T]) has 
its origin in the quantum fluctuation of the scalar in- 
flaton field 4>. The first-order perturbation S(f>, where 
4> = 4>o{t) + 5<fi(x), obeys the wave equation 



dfScf) + 3Hd t 5<j> - a 



- 2 E< 



> + m 2 (50 = O, (4) 



Let us focus first on these tensorial perturba- 
tions. The form of the modes is then h^(t, x) — 

(-167rGr7r/4(27r) 3 a 2 ) 1 / 2 i7i 1) (-fcr)e 4fe . The index of 
the Bessel function is v = \/9/4 + 3e and the so-called 
conformal time, r = J dt/a(t), is given by r = —(1 + 
e)/aH. Notice that, at early times, the amplitude of 
oscillations depends on k in a way similar to that of a 
massless field in Minkowski space. As time evolves, the 



comoving wavelength reaches the Hubble horizon length 
at t^. A few Hubble times after horizon exit the ampli- 
tude freezes completely and, for all subsequent times t, 
one has the constant value \h^\ 2 = GI ^. i ^ k " 1 . Because of 
the loss of phase information, the modes of the pertur- 
bations soon take on classical properties @. The freez- 
ing amplitude is usually codified through the quantity 
A 2 (k,t), defined in general by A 2 = 47rfc 3 |/i^(i)| 2 , and 
evaluated at the horizon crossing (or a few Hubble 
times after it). Taking into account the two polariza- 
tions, one finally obtains a nearly "scale free" tensorial 
power spectrum (J3J. In a similar way one obtains the 
scalar power spectrum ©. 

At this stage it is important to remark that the above 
definition of the power spectrum is such that the quan- 
tum fluctuations of the perturbations in position space 
satisfy the relation 



(h 2 (x,t))=l ?A 2 (fc,t) 



(5) 



This quantity represents the variance of the Gaussian 
probability distribution associated to h(x,t). However, 
in the form given by (|5|), it is divergent. It might be 
argued, as is common when dealing with random fields, 
that this divergence can be eliminated by smoothing out 
the field on a certain scale R to remove the Fourier modes 
with fc _1 < R. It can also be argued that the value of (h 2 ) 
is unimportant and regard the (finite) two-point function 
(h(x\)h{x2)) , uniquely defined by A^fc), as the basic ob- 
ject. However, it is our view to regard the variance as 
the basic physical object, which defines the amplitude of 
fluctuations in position space, and treat h(x, t) as a quan- 
tum entity. Therefore, if one regards h{x 1 1) as a quantum 
field, and {h 2 (x,t)) its two-point function at coincident 
points, the divergences must be handled in a different 
way. Taking into account the large k asymptotic form of 
the modes one can estimate the form of the integral at 
an arbitrary time 



where a(t) is the expansion factor of the unperturbed 
homogeneous and spatially flat metric ds 2 = —dt 2 + (h 2 (x,t)} 
a 2 (i)<£z; 2 . The effective mass term, which is necessarily 
small in the slow-roll approximation, is given by the sec- 
ond derivative of the potential: m 2 = V"(4>q). Moreover, 
the fundamental relation ([3]) has also the same quantum 
origin. The two independent polarizations of tensorial 
modes can be described by a couple of scalar fields /i+ jX 
obeying the above wave equation with zero mass. The 
relation between h (we omit the subindex + or x) and 
8<p is given by 84> = h/ VWttG. 



dk 16nGk 3 
k At: 2 a 3 



a , 



(2 + 3e) 
2fc 2 r 2 J 



(G) 



The first term produces a quadratic divergence, which is 
the typical singularity of a quantum field in Minkowski 
space. The second term produces a logarithmic diver- 
gence, which is typical of a massless field in an expand- 
ing universe. For a quantum mechanical system, with 
a finite number of degrees of freedom, such a divergent 
behavior would not arise and one need not worry about 
the definition of the physical power spectrum. However, 
a quantum field is neither a random field nor a quan- 
tum mechanical system with a finite number of degrees 
of freedom. The existence of divergences tells us that 
the nature of a quantum field is more involved. If we 
wish to have a finite expression for (h 2 (x,t)} one should 
subtract the divergences in a consistent way. In tech- 
nical words, one should renormalize the quantum field. 
Since the physically relevant quantity (power spectrum) 



3 



is expressed in momentum space, the natural renormal- 
ization scheme to apply is the so-called adiabatic sub- 
traction [Ioj |. as it renormalizes the theory in momen- 
tum space. Adiabatic renormalization 0, d, [H[ removes 
the divergences present in the formal expression ([5]) by 
subtracting (167rGfc 3 /47r 2 a 3 )[l/w + a 2 /2a 2 w 3 + a/2aw 3 } 
mode by mode in the integral (0, where w is the fre- 
quency of the modes. The subtraction of the first term 
(167rGfc 3 /47r 2 a 3 w) to cancel the typical flat space vacuum 
fluctuations was already considered in [l^]. However, the 
additional terms, proportional to a 2 and a are necessary 
to properly perform the renormalization in an expanding 
universe. Subtracting consistently the divergent terms 
one obtains the finite expression 



(h 2 (x,t)) 



dk l6nGk 3 , ttt 



o 

a r 



k At: 2 a 3 
, (2 + 3e) 

2k 2 T 2 1 



a^\H^\-krt 



(7) 



For the idealized case of a strictly constant H the sub- 
traction exactly cancels out the vacuum amplitude [To| 

since H^ 2 {x) = i exp(ix)^/2/7rx(— 1 +i/x). From this it 
follows that in pure de Sitter inflation there is no produc- 
tion of gravitational waves. One would thus expect that 
the tensorial power spectrum were proportional to the e 
slow-roll parameter, which parameterizes the slow change 
in H(t). The fluctuations still acquire classical properties 
through dccoherence in a time of order H~ x after their 
wavelengths become larger than the Hubble radius dur- 
ing inflation 0. Therefore, it is natural to evaluate the 
power spectrum (with the corresponding adiabatic coun- 
terterms) a few Hubble times after the time t^. Since 
the results will not be far different from those at tk, we 
use that time to characterize the results. We comment 
further on this point after (fT4|) . Evaluation of the above 
integral at the time tk leads to the redefinition 



A 2 (fc) = AGH 2 



1 



|tfW(-fcr = l + e) 



(4 



2tt 



(8) 

which turns out to be proportional to e, A^fc) = 
al6irG(H(tk)/2Tr) 2 e(tk) + 0(e 2 ), where a is a numerical 
coefficient of order unity, a w 0.904. Taking into account 
the usual conventions, we get, at leading order in e, 



P t (k) = — y , 
V ' M 2 V 



8a (H{t k ) 



2tt 



e(t k 



(9) 



One can proceed in the same way to reevaluate the 
scalar power spectrum. However, to be precise in the 
calculation, one must keep the mass term in (j4]) and 
also take into account the slow decay of the Hubble 
rate H, as above. The latter is controlled by the slow- 
roll parameter e, while the former is captured by the 
slowly changing parameter 77 = Mp(V" /V), which gives 
m 2 = 3ijH 2 . The form of the scalar modes is the same 
as that for the tensorial ones, up to the fact that the 
coefficient y/lQirG is now absent and the Bessel index is 



now v = 9/4 + 3e — 3r/. The adiabatic subtracting terms 
are then those already present for a massless field, with 
w(t) = \/k 2 1 a 2 + m 2 , together with new terms propor- 
tional to the mass (k 3 /47r 2 a 3 )[m 2 d 2 /2a 2 w 5 +m 2 a/4aw 5 — 
5m 4 d 2 /8a 2 w 7 ]. Therefore, after renormalization, the fi- 
nal result for the scalar power spectrum is 



H(t k 
2M 2 e(t k ) V 2ir 



(ac(tfe) + 3^(tjk)) , (10) 



where (3 « 0.448 is another numerical coefficient. In the 
case of eternal de Sitter expansion, one recovers the pre- 
vious results of [l(| EH • 

The above expression, together with © for the tenso- 
rial power spectrum, constitute our main results and lead 
to a change in the testable predictions of inflation. To see 
this we have to consider the scalar and tensorial spectral 
indices. The standard expression for them, in terms of 
the slow-roll parameters, is n s — 1 = — 6e + 2-q , n t = — 2e. 
However, if we invoke renormalization we get 



n t = 2(e - 1]) 



(11) 



-6e+2r]- 



4ae 2 + (6/3 - 2a)er] - 30(n' t /2 - 4e(e - rj)) 
ae + 3/3?7 



( 12 ) 

where n' t is the running of the tensorial index n' t = 
drit/d\nk. Note that n' t can be expressed in terms of 
the slow-roll parameters as n' t = 8e(e — 77) + 2£, where 
£ is defined by £ = M P (V'V" /V 2 ). The above formu- 
lae provide, implicitly, algebraic relations between phys- 
ical observables. The most important one is the relation 
between the tensor-to-scalar ratio r = P t /Bji with the 
spectral and running indices: r = r(n t ,n s ,n' t ). In the 
simplest case of n' t = 0, and taking the approximation 
a « 2/3, the new consistency condition becomes 



96 11 f1 ^ 96 . 
l-n s + -n t + Y] J(l-n s ) 2 + -n 2 



(13) 



This relation could be potentially checked, and compared 
with the standard one r = Sn t , in near future high- 
precision anisotropy measurements of the CMB. How- 
ever, we can already contrast, partially, our predictions 
with the standard ones on the basis of the five year 
WMAP results by performing a model by model anal- 
ysis. As a representative example we shall reexamine the 



(4-p), 



> p and corn- 



monomial chaotic potential V(tf>) = XM 
pare with standard theoretical and experimental results 
H- It is not difficult to obtain that 



ap 



(3j9(p-l)/2 + ap/4)JV ' 



1 - n.. 



4' 



where N = hxa en d/o,wMAP is the number of e- folds 
of inflation between the end of inflation and the epoch 
when the wavelength of fluctuations that WMAP detects 
left the horizon. If the adiabatic counterterms were 
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Figure 1: Plot of r versus n s . The contours show the 
68% and 95% CL derived from WMAP5 (in combination 
with BAO+SN) @]. We consider V{<j>) = m 2 <j> 2 (solid line), 
V((f>) = X(f> 4 (dashed line). The symbols show the prediction 
from each of this models in terms of the number N of e-folds. 
The top part corresponds to the prediction of our formulae, 
while the bottom one corresponds to the standard prediction. 



evaluated some e-folds n after tk but still well before 
the end of inflation (ne <C 1), then the value of a in 
the tensorial power spectrum would change to a ~ 2n. 
A similar calculation for the scalar power spectrum 
should also be carried out. This computation should 
deal with the gauge-invariant quantity IZk, which is 
conserved outside the Hubble sphere and satisfies the 



same equation as /ig, up to the replacement a — > aifio/H. 
The result is then a w 6n, 3/3 w — 2n. One then finds 
that the spectral index n s in (fH)) is unchanged, and the 
ratio r becomes r = Ap 2 /(p + 2)N, which is not sensitive 
to the value of n. If the counterterms are evaluated far 
beyond the end of inflation, one recovers the standard 
predictions r = Ap/N and 1 — n s = (p + 2)/2N . However, 
since primordial fluctuations acquire classical properties 
through decoherence when their wavelengths become 
larger than the Hubble radius @, we find it natural to 
evaluate the spectra at a time close to tk, which implies 
significant deviations from the standard prediction (see 
(fT4|)). Note that, with the present understanding of 
the nonlinear aspects of quantum gravity, it is difficult 
to reach a definitive conclusion regarding this question, 
so the fact that there are observable differences offers 
a deep way to experimentally probe this question. In 
fact, if we compare the results at tk with the WMAP 
5-year data for the representative values p = 2 and p = 4 
(Figure 1) we find that both models are compatible 
with the experimental data for the reasonable range of 
N between 50 and 60 (Figure 1 top). This is in sharp 
contrast with the prediction of the standard approach 
(Figure 1 bottom) where the monomial potential with 
p = 4 is excluded convincingly. With the new predictions 
both models are compatible with WMAP data for all 
the values of N between 50 and 60. This is also true 
even if the counterterms are evaluated n e-folds after 
tk, with ne < 1. The alternative predictions presented 
here may soon come within the range of measurement 
of high-precision CMB experiments, as expected in the 
PLANCK satellite mission. Acknowledgements 
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